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orbit  and  tracking  data  to  determine  errors  in  station  loca¬ 
tion.  It  also  documents  a  Philco  2000  Computer  program 
which,  given  the  standard  deviations  in  the  tracking  data, 
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FOREWORD 

Technical  Note  WDL-TN62-14  has  been  prepared  by  the  Philco  WDL 
Mathematical  Analysis  Department  for  submittal  to  AFSSD  for  information 
purposes.  This  Technical  Note  is  within  the  scope  defined  by  Paragraph 
4.2.1,  AFBM  Exhibit  58-1,  "Contractor  Reports  Exhibit,"  dated  1  October 
1959,  as  revised  and  amended. 

The  material  presented  in  the  Technical  Note  was  developed  in 
conjunction  with  Tracking  Simulation  and  Evaluation  and  Advanced 
Trajectory  Analysis  Studies  conducted  by  Philco  WDL  under  Exhibit  "A" 
of  Definitive  Contract  AF04(695) -113 ,  and  Paragraph  1.2. 1.2  of  AFSSD 
Exhibit  61-27A,  "Satellite  Control  Subsystem  Work  Statement,"  dated 
15  February  1962. 
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SECTION  1 
INTRODUCTION 


Current  and  projected  mission  requirements  place  a  severe  demand 
on  the  ability  to  accurately  determine  the  position  of  a  satellite. 

Since  the  equations  of  motion  for  a  satellite  are  formulated  in  inertial 
space  and  the  tracking  information  is  given  in  local  coordinates,  it  is 
necessary  to  know  the  location  of  the  tracking  stations  to  an  accuracy 
compatible  with  the  inherent  noise  level  of  the  tracking  equipment. 

The  existence  of  tracking  data,  and  hence  an  orbit  for  a  satellite, 
can  be  used  to  determine  the  location  of  a  station  in  relation  to  other 
stations. 

The  observational  data  available  from  a  satellite  which  is  used 
for  the  purpose  of  locating  a  station  on  the  geoid  are  the  elevation 
(E),  azimuth  (A) ,  slant  range  (S),  and  range  rate  (S).  Not  all  of 
these  quantities  will  necessarily  be  available. 

A  camera  at  the  site  would  give  angular  data  about  the  orbit;  an 
antenna  might  give  range  data  in  addition  to  giving  angular  data.  A 
computer  program  has  been  written  which  indicates  the  errors  in  station 
location  resulting  from  errors  in  the  camera  or  antenna  data.  It  is 
assumed  the  errors  in  azimuth  and  elevation  (and  range  as  a  possible 
option)  are  randomly  distributed  with  zero  mean  and  specific  variances. 
As  a  prediction  of  a  survey's  accuracy,  the  expected  variances  in  the 
data  can  be  used  to  find  the  resulting  covariance  matrix  in  station 
coordinates.  The  eigenvalues  of  the  matrix  are  also  calculated. 
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EQUATIONS 


An  elliptical  orbit  is  assumed  in  the  computations.  The  orbit  is 
defined  by  the  elliptic  elements  tQ,  o>,  e,  a,  i,  and  ft. 


Over  each  pass  q  is  obtained  by  an  iterative  process  from' 


t(n)  “  t(nQ) 


2  tan 


-1 


1  + 


e  sin  A 

7^ 


e  +  e  cos  A 


ey  1-e2  sin  A 
1  +  e  cos  A 


A  -  T]  -  ojq 

A  -  Tj  -  0) 

'o  o 


In  this  equation 

t(q)  -  time 

-  eccentricity 

-  n2a3  -  62625.53 
1 

a(l-e2) 

These  elliptic  elements  are  then  transformed  to  inertial  cartesian 

coordinates  and  thence  to  local  polar  coordinates  by  standard  transforma- 
2 

tlon  equations  .  The  local  polar  coordinates  are  then  used  to  compute 
the  elements  of  the  covariance  matrix. 


e 

L 


1.  E.  D.  Callender,  "A  Solution  of  the  Equations  of  Motion  for  a 
Near-Earth  Satellite,"  to  be  published. 

2.  "Program  I  Simulation  and  Evaluation  Study  Report,".  Philco 
Corporation,  WDL-TR1596,  1  October  1961. 
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If  the  station  coordinates  are  found  by  a  differential  corrections 
technique  using  least  squares,  the  covariance  matrix  will  be  given  by 


where  A  is  the  weighted  Jacobian  matrix  and  the  summation  is  over  all 

3 

times  at  which  data  is  provided  .  The  use  of  A  for  both  azimuth  and 
for  the  Jacobian  matrix  is  standard  notation  and  should  cause  no  con¬ 
fusion.  The  elements  of  A  are  as  follows  (notation  is  0q  ■  longitude,. 
0  -  latitude,  h  -  height  above  the  geold,  and  Rg  -  radius  to  station): 

cos  9  (x  sin  0  -  y  cos  0) 

(x  sin  9  cos  0  +  y  sin  0  sin  0  -  z  cos  9) 

(Rg  -  x  cos  9  cos  0  -  y  cos  9  sin  0  -  z  sin  9). 


90  S 


as  .  I 

3h  S 


£s 

80 


£S 

8h 


~  f  +  IT  COB  6  {f*  8in  0  "  *  cos 


9s  S 
89  S 


+  .0  (y  sin  0  +  x  cos  0)J 

^  ^sin  0  g,  cos  0  +  y  sin  0) 

+  0  (y  cos  jT  -  x  sin  0)J  -  z  cos  0^ 


!§  i.l  (cos  0  (x  cos  0  +  y  sin  0 ) 

8h  s  s  L 

+  z  sin  0  +  0  cos  0  (y  cos  0  -  x  sin  0)J 


3.  "A  User's  Manual  for  Three  Tracking  Simulators,"  Fhilco  Corporation,. 
WDL-TN62-1. 
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c. 


£e  _  1 

c)0  S  cos  E 


0 


(y  cos  0 


x  sin  0)  cos  0 


8S 

90 


sin  E 


0 


3E  1 

—  m  «— — «— i — 

3q  S  cos  E 


|je  cos  9  -  (x  cos  0  +  y  sin  0)  sin  9 


ds 

de 


sin  E 


0 


9E 

9h 


1 

S 


£s 

8h 


tan  E 


1  1 
S  cos  E 


d; 


A 


tan 


-1 


N 

D 


So,  in  general. 


8^  9d 
JP..9fl.  r..JLja. 

t>2 


N  -  -  x  8 in  0  +  y  cos  0 
D  -  -  sin  0  (x  cos  0  +  y  sin  0)  +  z  cos  0 


(1) 


(2) 


(3) 


£n 

90 


x  cos  0  -  y  sin  0 


—  -  -  sin  0  (y  cos  0  -  x  sin  0) 
00 


9n 

9e 

do 

Be 

8a 

8h 


0 

-  cos  9  (x  cos  0  +  y  sin  0)  -  z  sin  9 
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2-~3 


PHILCO 


WKITHN  DEVELOPMENT  LABORATORIES 


WDL-TN62-14 


In  the  preceding  equations,  0  -  flt  +  0  q  +  A  0 ,  where  0  is  the  rate 
of  the  earth's  diurnal  rotation,  0q  is  the  longitude  from  Greenwich, 
and  A  0  is  the  phase  angle.  However,  we  have 

8,  _8s 

h  a,  a, 

O  O 

so  that  the  preceding  expressions  can  be  used  directly  in  the  covariance 
matrix. 
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SECTION  3 
COMPUTER  RESULTS 

Two  orbits  were  used  as  test  cases  for  the  program:  6161  and  6l£t81 
the  first  having  an  eccentricity  of  approximately  0.1  and  a  perigee  of 
469  statute  miles  and  the  second  having  an  eccentricity  of  0.01  and  a 
perigee  of  2180  statute  miles.  Palo  Alto  was  used  as  the  observation 
station. 

For  the  6151  orbit,  three  passes  were  observed,  two  overhead  and 
one  at  a  maximum  elevation  of  approximately  4-1° .  The  time  interval  for 
the  two  overhead  passes  were  approximately  30  minutes  duration,  while 
the  pass  at  the  lower  elevation  angle  was  about  24  minutes  in  length. 
With  a  10-second  interval  between  observations,  this  allowed  about  180 
observations  for  the  30-minute  pass  and  about  144  observations  for  the 
24-minute  pass. 

These  passes  were  used  in  combinations  of  two  to  obtain  four  sets 
of  results  using  a  standard  error  for  azimuth  and  elevation  of  5  milli- 
radlans  for  3  sets  and  3  milliradians  for  the  remaining  set. 

For  the  610(81  orbit,  two  passes  were  used  both  singly  and  in  com¬ 
binations  with  equal  errors  in  elevation  and  azimuth  of  2  milliradians 
and  5  milliradians,  thus  obtaining  6  sets  of  results  for  the  covariance 
matrix. 

The  listing  of  the  orbit  and  computer  results  is  given  in  Appendix 
A,  and  a  description  of  the  Computer  Program  and  ALTAC  Listing  is  given 
in  Appendix  B. 
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ORBIT  AND  COMPUTER  RESULTS 
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DESCRIPTION  OF  COMPUTER  PROGRAM  AND  ALTAC  LISTING 
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INPUT  FOR  DC11  PROGRAM 
(GEODETIC  SURVEY) 

One  complete  set  of  data  consists  of  at  least  6  data  cards.  Each 
card  Is  divided  into  8  fields,  10  slots  long. 


SUMMARY  OF  DATA: 


CARD  1 


CARD  2 


CARD  3 


CARD  4 


TO 

t 

0 

(MIN)  > 

ATAO 

^0 

(DEG) 

INITIAL  CONDITIONS 

OHL 

(0 

(DEG) 

OF  THE  SATELLITE; 

ECC 

“ 

>  (FLOATING  POINT  NUMBER, 

AX 

a 

(N.M. ) 

x-xxx+x) 

FI 

1 

(DEG) 

OHB 

a 

(DEG) 

J 

NOOBPASS 

NO.  OF  OBSERVATIONAL  PASSES 

(A  FIXED  POINT  NUMBER) 

FLAT 

9 

(DEG)  " 

FLONO 

*E 

(DEG  E) 

)  STATION  COORDINATES 

ALT 

h 

(FT) 

DELPHIO 

A  0 

(DEG) 

DIMEN 

DIMENSION  OF  COVARIANCE  MATRIX 

(2  -  LAT  &  LONG,  3  -  LAT,  LONG, 

&  ALT) 

(A  FIXED  POINT  NUMBER) 
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a 


CARD  5 


CARD  6 


SIGS 

SIGD 

JFLAG1 

SIGA 

SIGE 

T1 

T2 

DELT 


U 

s  (N.M. ) 

0  (FT/SEC) 


VSEC)  'j  (MI. RAD) 

A  l  IF  \  IF 

JFLAG1  -  OJ  JFMG1  “ 


tl  (MIN) 
t2  (MIN) 
At  (SEC) 


I  REPEAT 
1  /  FOR 
EACH 

OBSERVATIONAL 

PASS 


J 


I  I  -  0,  NO  MORE  DATA  CASES 

I  -  I,  ANOTHER  COMPLETE  SET  OF  DATA 

FOLLOWING  STARTING  WITH  CARD  1 


T 


1 


i 

i 
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E 
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IS 


[ 

1" 

I 

I 

I 

I 


c 
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SAMPLE  SET  OF  DATA 
FOR  DC11 


10 

20 

30 

40 

50 

60 

70 

0.44+3 

0.0+1 

0.27+3 

0.1031-0 

4.3+3 

0.4+2 

0.19+3 

1 

0.37+2 

2  .-37+2 

0.201-2 

0.16+2 

2 

0.0+1 

0.0+1 

0 

0.1+2 

0.10+2 

0.45+3 

0.48+3 

0.10+2 

1 

>CA: 


SE  1 


0.44+3 

0.0+1 

0.27+3 

0.1+0 

4.3+3 

0.4+2 

0.2+3 

1  1 

2 

. 

1  1 

0.37+2 

2.37+2 

0.2+2 

0.16+2 

1  1 

3 

0.1+1 

0.2+1 

1 

0.01+1 

0.02+1 

0.4+3 

0.5+3 

0.1+2 

1 - 

0.0&+1 

0.1+1 

1 

0.01+1 

0.17+1 

0.4+3 

0.2+2 

-  ~ _ 

0 

CASE  2, 

2  PASSES 
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VARIABLE 

NAME  &  DESCRIPTION 

t 

o 

TO 

TIME  OF  INITIAL  CONDITIONS  (MIN) 

no 

ATAO 

V  +  CD  (DEG) 

(0 

OHL 

ANGLE  OF  PERIGEE  (DEG) 

€ 

ECC 

ECCENTRICITY 

a 

AX 

SEMI-MAJOR  AXIS  (N.M. ) 

1 

FI 

ANGLE  OF  INCLINATION  (DEG) 

n 

OHB 

OMEGA  (DEG) 

e 

FLAT 

LATITUDE  (DEG) 

0 

FLONO 

LONGITUDE  (DEG) 

h 

ALT 

ALTITUDE  (FT) 

60 

DELPHIO 

ANGLE  BETWEEN  INITIAL  X-AXIS 

AND 

O 

GREENWICH  AT  T  ZERO 

ffSR 

SIGS 

<J  "N 

SLANT  RANGE  (N.M.) 

aSD 

SIGSD 

°DOPPLER  (FT/SEC) 

STANDARD 
DEVIATIONS 
i  OF  NOISE 

or 

A 

SIGA 

^AZIMUTH  (MILLIRAD  OR  SEC) 

IN  THE 
TRACKING 

or 

E 

SIGE 

^ELEVATION  (MILLIRAD  OR  SEC) 

DATA 

J 

h 

T1 

TIME  IN  (MIN) 

t2 

T2 

TIME  OUT  (MIN) 

At 

DELT 

TIME  STEP  (SEC) 
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COMPUTED  CONSTANTS 


C200 

90.0 

C100 

2.0/a 

C2 

a/180 

CIO 

180/a 

Cll 

o 

H 

o 

to 

C101 

1.0/C10 

C3 

(180) 3600 /a 

Clio 

1.0/C100 
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ATAO 

SUMMARY  OF  EQUATIONS 

-  C2  *  ATAO 

OHL 

m 

C2  *  OHL 

V 

m 

ECC  *  cos  (OHL) 

WA 

as 

ECC  *  sin  (OHL) 

TO 

■i 

60.0  *  TO 

FL 

■ 

1.0/ (AX  *  (1.0  -  ECC 

ALT 

m 

ALT/6076.10333 

FLAT 

m 

C2  *  FLAT 

FLONO 

m 

C2  *  FLONO 

DELPHIO 

m 

C2  *  DELPHIO 

COS.l 

m 

cos  (FLAT) 

SIN  1 

- 

sin  (FLAT) 

DELPHIQ 

- 

.  C100  *  DELPHIO 

FMEANMO 


62625.53 
,3 


AX 


C4 


- 


ECC 


FATAO  -  - 


2-0  *  tan"1  *  SrlnJilr°:0WL) 

\  1.0  +  C4  +  ECC  *  cos 


(ATAO-OHL) t 


ECC  *  C4  *  sin  (ATAO-OHL) 
1.0  +  ECC  *  cos  (ATAO-OHL) 


FLATA  -  C100  *  FLAT 

FLONA  -  C100  *  (FLONO  -  .00007292*TO> 

T1  -  60.0  *  T1 

T2  -  60.0  *  T2 
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E 

A 

FI 

OHB 

COSE 

SINE 

COSA 

SINA 

FLON 

COS  2 

SIN2 

SLON 

SLAT 

SALT 

SDLON 


SDLAT 


S  DALT 


SUMMARY  OF  EQUATIONS  (CONTINUED) 

CllO  *  E 
CllO  *  A 
C200  *  FI 
C200  *  OHB 
cos  (E) 
sin  (E) 
cos  (A) 
sin  (A) 

FLONO  +  DELPHIO  +  O.OOOQ7292  *  (T  -  TO)  _ _ _ 

cos  (FLON) 
sin  (FLON) 

R  *  C0S1  *  (X  *  SIN2  -  Y  *  C0S2)/S 
R  *  (X  *  SIN1  *  C0S2  +  Y  *  SIN1  *  SIN2  -  Z  *  C0S1)/S 
(R-X  *  C0S1  *  COS 2  -  Y  *  C0S1  *  SIN2  -  Z  *  SIN1)/S 

-  SLON  *  SD  +  R  *  C0S1  * 

S 

[XD  *  SIN2  -  YD  *  C0S2  +  FLOND  *  (Y  *  SIN2  +  X  *  C0S2)J  /S 

-  SLAT  *  SD  +  R  *  [siNl  * 

S  L 

^XD  *  COS 2  -  YD  *  SIN2  +  FLOND  *  (Y  *  C0S2  -  X  *  SIN2)} 
-ZD  *  COSlj  /S 

-  SDALT  *  SD  -  [XD  *  COS2  *  C0S1 

S  «- 

+  YD  *  SIN2  *  C0S1  +  ZD  *  SINl  +  FLOND  *  COS1  * 

(Y  *  COS 2  -  X  *  SIN2)]  /S 
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ELON 

TANK 

ELAT 

EALT 

Cl 

ALON 

ALAT 

AALT 


SUMMARY  OF  EQUATIONS  (CONTINUE I;) 


[COSI  *  (Y  *  COS  2  -  X  *  SIN2)  -  SINE  *  SI.On]/ 


(S  *  COSE) 


-  SINE /COSE 

-  Qz  +  COSI  -  SIN1  *  (X  *  COS 2  +  Y  *  SIN2) 

-  SLAT  *  SINE]  /  (S  *  COSE) 

-  -  SALT  *  TANE  -  LjJ)  *  COSE 

S  S 

-  1.0/  (FN2  +  D2) 

-  Cl  *  [D  *  (-X  *  COS2  -  Y  *  SIN2)  +  FN  *  SINl  * 

(Y  *  COS 2  -  X  *  SIN2)] 

-  Cl  *  FN  *  [COSI  *  (X  *  COS 2  +  Y  *  SIN2)  -  Z  *  SIN2] 

■  .  0.0 


P  MATRIX 

SLAT 

SLON 

SALT 

SDLAT 

SDLON 

SDALT 

ALAT 

ALON 

AALT 

ELAT 

ELON 

EALT 
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SUMMARY  OF  EQUATIONS  (CONTINUED) 
FOR  2  DIMENSIONAL  PROBLEM: 

+  V1  P2l  P21 

+  W2  P21  P22 

+  W22  P22  P2I 
.  2  „  „ 


COMAT n 

m  M  ^  P  ! 
1  *11  ‘ 

COMAT 12 

m  wl2pll: 

comat21 

-  Hl2p12; 

comat22 

2 

at  TJ  P 

w  1  *12 

C5  - 

COMAT n  COMAT, 

C5  - 

1.0/C5 

®U  “ 

C5  *  COMAT22 

<*12  " 

-C5  *  C01AT21 

“21- 

-C5  *  COMAT 12 

CM22  " 

C5  *  COMAT11 

22 


12 


21 


c5i  -  *  (4*°*c5> 

EIGENI  -  VK“u + ®22) + csii  *  °-5 


EIGEN 2  -  \/[(cMll  +  «2*)  '  C51]  *  °*5 


CONVERT  CM  ARRAY  TQ  DEGREES  SQUARED: 


CMU 

-  Cll 

<CMn) 

«12 

-  Cll 

(cm12) 

cm2i 

-  Cll 

(cm21) 

0*22 

-  Cll 

C«22> 
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SUMMARY  OF  EQUATIONS  (CONTINUED) 
FOR  3  DIMENSIONAL  PROBLEM: 


COMAT .. 

m 

»,2 

P11 

pll  + 

w  2 
W2 

P2I 

P21  + 

w  2 
W3 

P31 

P31 

+  W42 

P41 

P41 

COMAT12 

- 

.•l2 

P11 

PI2  + 

w  2 
W2 

P21 

P22  + 

w  2 
W3 

P31 

P32 

+  W42 

P41 

P42 

COMAT 13 

SI 

“l2 

P11 

P13  + 

V  2 
W2 

P21 

P23  + 

w  2 
3 

P31 

P33 

+  W42 

P41 

P43 

COMAT21 

B 

P12 

pll  + 

w  2 
W2 

P22 

P21  + 

w  2 
W3 

P32 

P31 

+  W42 

P42 

P41 

COMAT22 

“,2 

P12 

P12  + 

w  2 
2 

P22 

P22  + 

w  2 
W3 

P32 

P32 

+  W42 

P42 

P42 

COMAT23 

Si 

V  2 
W1 

P12 

P13  + 

w  2 
"2 

P22 

P23  + 

w  2 
W3 

P32 

P33 

+  W42 

P42 

P43 

C(»1AT31 

SI 

w  2 
W1 

P13 

Pll  + 

w  2 
W2 

P23 

P21  + 

«  2 
W3 

P33 

P31 

+  W42 

P43 

P41 

COMAT32 

m 

V  2 
W1 

P13 

P12  + 

w  2 
"2 

P 

r  23 

P22  + 

H  2 
W3 

P33 

P32 

+  W42 

P43 

P42 

COMAT33 

- 

w  2 
W1 

P13 

P!3  + 

H  2 
W2 

P23 

P23  + 

w  2 
W3 

P33 

P33 

+  W42 

P43 

P43 

C20  - 

COMATn 

*COMAT22*COKAT. 

33  + 

2.0  * 

C0MAT12' 

*comat23*comat 

31 

-comat132*comat22  -  comat232*comat11 
-comat122*comat33 


C21 

m 

1.0/C20 

CMH 

m 

c21*(comat22 

COMAT33  -  COMAT23  ) 

CMi2 

m 

-C21*(COMAT21 

COMAT33  -  COMAT23  C0MAT31) 

cm2i 

m 

CMi2 

CMi3 

m 

C21*(COMAT21 

COMAT32  -  C0MAT31  COMAT22) 

cm3i 

m 

CMi3 

CM22 

m 

C21*(  COMAT  1]L 

COMAT33  -  COMAT132) 
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SUMMARY  OF  EQUATIONS  (CONTINUED) 


e 


FOR  3  DIMENSIONAL  PROBLEM  (CONT'D.): 


CM33 

cs 

C21*(C0MAT11  C0MAT22  -  COMAT^ 

“23 

B 

-C21*(C0MA,T11  COMAT32  -  COMAT12 

“32 

“ 

“23 

GOING 

TO 

SUBROUTINE  FOR  EIGEN  VECTOR  AND 

CMu 

- 

Cll  (CMn) 

“12 

- 

cn  <cm12) 

CMi3 

- 

cio  (cm13) 

“21 

B 

cn  (cm21) 

“22 

m 

cu  <cm22) 

“23 

- 

cio  (cm23) 

“31 

a 

CIO  (CM31) 

“32 

- 

CIO  <cm32) 

“33 

“33 

C5  ' 

- 

COMAT X1  COMAT22  -  COMAT x 2  COMAT, 

C5  ■ 

- 

1.0/C5 

“ll 

B 

CS  *  COMAT22 

“12 

B 

-C5  *  COMAT21 

“21 

- 

-C5  *  COMAT 12 

“22 

B 

C5  *  COMATn 

21 
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n 

SUMMARY  OF  EQUATIONS  (CONTINUED)  j'j 

AFTER  GOING  TO  SUBROUTINE  FOR  EIGEN  VECTOR  AND  MATRIX  (CONT.'D.):  n 


cm21  -  cii  (cm21) 
cm22  .  CU  (CM22) 


1 

B 

0 

0 

0 

11 

0 

D 
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407  LISTING  (CONT'D) 
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DISTRIBUTION  LIST 


Address  Number 

Commander 

Space  Systems  Division 
Air  Force  Systems  Command 
Air  Force  Unit  Post  Office 
Los  Angeles,  California 

USAF  Contract  Support  Detachment  No.  3 
Philco  Corporation 
Western  Development  Laboratories 
Palo  Alto,  California 

ASTIA 

Arlington  Hall  Station 
Arlington  12,  Virginia 

Philco  Corporation 

Western  Development  Laboratories 

Palo  Alto,  California 
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